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MTIOML ADVISORY COMMITTEE FOR AEROMUTICS 
RESEARCH MEMORANDUM 

GENERAL THEORY OF WAVE-DRAG REDUCTION FOR COMBINATIONS 
EMPLOYING QUASI -CYLINDRICAL BODIES WITH AN APPLI- 
CATION TO SWEPT-WING AND BODY COMBINATIONS 
By Jack N. Nielsen 

SUMMARY 


The wlng-hody interference theory of NACA TN 26 tT applied to symmet- 
rical wings in comhination with quasi -cylindrical bodies permits the 
direct calculation of pres sure -distribution changes produced by body 
shape changes. This theory is used to. determine the relative magnitude 
of the wave-drag reduction produced by changes in cylinder cross-sectional 
area and that produced by changes in cross-sectional shape (without change 
in area). The body distortion is expressed as a Fourier series, and an 
integral equation is derived for the body shape for minimtim drag for each 
Fourier component. Thus the wave-drag reductions for the various Fourier 
harmonics are independent and additive. 

To demonstrate the use of the method, it is applied to. a body with 
a 45 ° swept wing having a chord of three body radii, a span of six body 
radii, sonic leading and trailing edges, and a biconvex airfoil section. 
The axially symmetric harmonic of body distortion caused a reduction of ■ 
wave drag of 22 percent of the drag of the wing panels on the cylindrical 
undistorted body. However, the second harmonic, which causes no change 
in volume, produced a drag saving of 42 percent. The wave -drag saving 
was obtained by decreasing the wing-panel drag by twice as much as the 
body drag was increased. The wing panels have negative wave drag for two 
harmonics. These reductions are accompanied by pressure gradients that 
are adverse for the boundary layer. With regard to fuselage shape, volume 
is removed from the sides of the fuselage and added to the top and bottom. 
A net volume less than the wing volume is removed if the distortions are 
limited to the interval of the root chord. Althou^ the second hanaonic 
removes no volume, it causes indentations about five times as deep as 
those due to the axially symmetric harmonic. As a result the wing thick- 
ness for -vdiich the full theoretical distortions can be utilized is limited 
to about 5 percent. 

A very simple theory for the axially symmetric harmonic can be 
derived by neglecting three-dimensional effects. This theory is valid 
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for small chord lengths. However, an approximate treatment of the three- 
dimensional effects yields a fairly simple theory applicable to longer 
chord lengths . The effects of semi span -radius ratio and chord-radius 
ratio on indentation shape, indentation volume, and drag saving have been 
determined using the approximate theories for the axially symmetric 
harmonic . 


INTRODUCTION 


The experimental verification by Whitcomb (ref. l) of the area rule 
has opened up a new realm of possibilities of reducing the wave drag of 
complete configurations at transonic and supersonic speeds. The theo- 
retical basis for the area rule lay undiscovered for some time in the 
Hayes drag formula (ref. 2), which contains the implication that by 
indenting a body, a saving in wave drag equal to the drag of the wing 
panels can be produced. Thus, in principle, the wave drag of a slender 
combination utilizing an infinite cylindrical body can be reduced to 
zero. More recently Lomax and Heaslet (ref. 3) have proposed the follow- 
ing method for determining wing -body combinations of minimum wave drag. 
Starting initially with a wing alone, they 

(1) use the invariance principle of Hayes to find distributions of 

multipoles along the longitudinal axis that represent the 
wing-alone flow field at large radial distances 

(2) designate the negatives of the multipole distributions as can- 

celing multipoles, place them along the longitudinal axis, 
and calculate a stream tube in the vicinity of the wing. 

The wing -body combination given by (2) has zero wave drag if all multi - 
poles are used. If, to obtain practical shapes, only a finite number of 
multipoles are used, the wave drag is a minimum for this ntmiber. 

In the foregoing method the problems of actual body shape and of 
pressure distributions are secondary considerations in the mathematical 
solution of the minimum -wave -drag problem. A method for minimizing wave 
drag of wing-body combinations dealing directly with the body shape and 
the pressure distributions has been developed based on the interference 
theoiy of references k, 5, and 6. The method, presented herein, is 
applicable to combinations with bodies which are approximately cylindri- 
cal in the region of the wing. There are no restrictions on the Mach 
number as long as it is greater than unity but not hypersonic. The 
principal restriction on the wing other than small thickness and sharp 
leading edges is that it possess a horizontal plane of symmetry or that 
its upper and lower surfaces be independent. It is the purpose of this 
paper to present in the first part a general theory of drag reduction 
for wings on quasi -cylindrical bodies. In the second part the method is 
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applied to the reduction in vave drag of a combination employing a swept 
wing with sonic edges. 


SYMBOLS 


All 


a 

Ag 

c 

Cr 

Cp 


^Dc 


^(C) 

®B(C)2n 

Dc 

I>S 

I>W(d) 

%(i) 

(i) 

%(i+d) 


(2) 

%(i+d) 


%(C) 

%(W) 

dx 


lengths are measured in body radii, 
body radius 

aspect ratio of exposed wing panels Joined together 
local wing chord or chord of rectangular* wing 
chord at wing-body Juncture 

wave -drag coefficient of combination based on exposed wing 
area , 

reduction in drag coefficient based on exposed wing area as 
a result of body distortion 

drag of body in combination 

drag of body due to 2nth Fourier harmonic 

drag of combination 

saving in drag of combination by body distortion 

drag of wing in combination with body due to body distortion 

drag of wing in combination with body due to interference 


drag of wing in combination due to interference and distor- 
tion pressures propagated along characteristic body cones 


drag of wing in combination due to interference and distor- 
tion pressures falling behind chai*acteristic body cones 

drag of wing panels on distorted body 

drag of wing panels acting as part of wing alone 

local slope of wing s'urface in streamwise direction 
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^‘2n(x) 

S2n(^) 

(i) 

So ^ 

&o ^ • 
h2n(x) 

^Oi(x), 

liOgCx), 

I 

m 

M 

n 

2n 
P ■ 

Pgn 

Poo 

p 

^B(d) 

PB(i) 

Pb(w) 

Pw(d) 


velocity amplitude function of 2nth hamonic neglecting tip 
effects 

• velocity amplitude function for a pair of opposing sources 
inclined at the Mach angle 

amplitude of body distortion due to 2nth harmonic 

, successive approximations to go(x) 

nonhomogeneous part of integral equation for ganCx) 
components of ho(x) 

upstream limit of body distortion 
slope of wing leading edge 
free -stream tfech ntanber 
normal to quasi -cylindrical body 
index of Fourier harmonics 
local static pressure 

half -amplitude of the 2nth Fourier component of the pres- 
sure acting on the body due to the wing alone neglecting 
tip effects 

free-stream static pressure 

p - p 

pressure coefficient, . — — 22 

body pressure coefficient due to body distortions 
body pressure coefficient due to interference 
body pressTire coefficient due to wing alone 
wing pressure coefficient due to body distortions 
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^w(i) 

%(w) 

q 

r,0 

R.P. 

S 

Si(x) 

So(x) 

Sp 

Sy 

u,v,w 

Vn 

vr 

V 

Vi 

Vw 

W2n(x,r) 

x,y,z 


Xl 


X2 


wing pressure coefficient due to interference 
wing pressvire coefficient due to wing alone 

1 p 

free -stream dynamic pressure, g pV 

polar coordinates in y,z plane; y = r cos 0, z = r sin 0 

radial distance to point on body of combination 

real part of a complex function 

semispan of wing-body combination 

local semispan of inner portion of wing plan form 
(sketch (c)) 

local semispan of outer portion of wing plan form 
(sketch (c)) 

area of one exposed wing panel 
area of two exposed wing panels 
maximum thickness of wing section 

perturbation velocities along x, y, and z axes, respectively 

velocity normal to quasi -cylindrical body, positive outward 

velocity along radius vector, positive outward 

free-stream velocity 

volume of body indentation 

volume of exposed wing 

influence function • 

coordinate axes with origin at vertex of wing alone, 

X measured downstream, y laterally to the pilot’s right, 
and z vertically upward 

intersection with body of most forward body characteristic 
touching exposed wing 

intersection with body of most rearward body characteristic 
touching exposed wing. 
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x»,y*,z» 


oblique coordinates; x' = x - my, y’ = y - mx, z’= z'^l - 


Y(x) 


Wo(l)d| 


Z(x) 


J |Wo(|)d| 


Zu 


'n 


'Pc 

CPi 

'Pw 

tip. 


local thickness of lower wing surface measured from chord 
plane, z = 0 

local thickness of upper wing surface measured from chord 
plane, z = 0 

- 1 

s 


^2n(x ^r)dr 


semiapex angle of wing section 

function defining the plan-form distribution of wing thick- 
ness (eq. (2)) 

constant^ Tvo = = 2 if n > 0 

free -stream density 

potential of distorted, body alone 

potential of complete combination 
interference potential, ~ 

potential of wing alone 

An asterisk (*) distinguishes a quantity associated with the wing 


I - GENEEIAL THEORY 

Statement of Problem and Determination of Boundary Conditions 

In references 4, 5^ S'^id 6 the boundary -value problem of determining 
the pressure distribution acting on a combination composed of a wing and 
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an infinite, circular, cylindrical. tody was solved. In the method, the 
flow field of the wing alone is pictTired as causing a virtual distortion 
of the body surface, and the pressiire field due to such distortion is 
computed term hy term for each Fourier heirmonic of the distortion. The 
pressure field on the tody is computed with the help of a table of char- 
acteristic functions or influence coefficients -vrtiich give the pressure 
coefficient at any point of the body due to a unit bump at the origin. 

By summing the pressure contributions of each bump on the body by numeri- 
cal integration, the entire pressure field is built up. An essentially 
different and more complicated method using axial distributions of 
sovirces and multipoles was used to obtain the pressure field off the body. 
However, it has recently been possible to extend the method for computing 
the pressures for points on the body to points off the body by expanding 
the tables of influence coefficients. The theory of the new. method is 
set forth in Appendix A. This method which has been applied to combina- 
tions with cylindrical bodies having virtual distortions is equally 
applicable to combinations having bodies with real distortions. Thus we 
can evaluate the pressure field and the presstire drag of a combination 
with a body of arbitrary small distortion using the above methods. The 
pressure drag of the combination can be minimized with respect to body 
shape by the calculus of variations to determine the best shape for mini- 
mum drag. In this part of the report the theory of determining this 
shape will be considered. 


In this paper the wing alone is taken to be the wing panels plus the 
portion of the wing blanketed by the body when the leading- and trailing- 
edge intersections with the body are joined across the body in- any arbi- 
trary manner. The wing -body combinations to which the analysis applies 
possess a vertical plane of symmetry and either a horizontal plane of 
symmetry or upper and lower wing surfaces that are indep>endent (e.g., a 
triangular wing with supersonic leading edges) . An additional assumption 
is that the body shape at any point does not greatly deviate from a cir- 
cular cylinder; that is, the body is quasi-cylindrical. Sketch (a) 
illustrates the winged portion of a 
wing-body combination together with the 
system of axes used in the analysis. 

All distances are based on a radius of 
unity so that the formulas in the anal- 
ysis can be written in simple form. 

The body cross section at any x posi- 
tion is taken as a circle on which are 
superimposed distortions that vary as 
cosines of even multiples of 9 to 
preserve a vertical plane of symmetry. 

Thus 


yj 



n=o 


( 1 ) 


Sketch (a) 
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The function g 2 n(x) is the shape amplitude function of the 2nth Fourier 
harmonic and is to he so determined that the wave drag is a minimum. 

With regard to the wing, no limitations are made on its plan form. The 
coordinates of the wing section at any spanwise station y are taken to 
he 


zu = -zz = tmTi(x - 3r + p,r) 


( 2 ) 


where is a function giving the thickness distribution and is zero off 
the wing. 

Consider now the boundary -value problem of determining the combina- 
tion potential. This problem is solved in essentially the same manner as 
in reference 4. For a cylindrical body the combination potential is the 
sum of the wing-alone and interference potentials, but when the body is 
distorted, an extra potential cpg must be added. 

cpc = + Ti + Tb (3) 


All potentials are to satisfy the vave equation in cylindrical coordinates. 

^ Sr ^ Sx^ 

The first houndaiy condition is that the derivative of the combination 
potential normal to both wing and body be zero with the body taken as 
the r = 1 cylinder and the wing as the z = 0 plane . Thus 

Scpp r = l O<0<2it (body) ^ 

-^ = 0 (4a) 

on d = 0,jt (wing) 

The second boundary condition is that there be no upstream waves. 


cpQ = 0 X < 0 


(4b) 


The second, boundary condition is satisfied naturally by using Laplace 
transfonn methods, but the first one is not^ so easily met. Consider 
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sketch (h) vhich shows a region of 
the quasi -cylindrical surface 
R(x,r,0) = constant in the neigh- 
borhood of point P, The direction 
cosines of the normal to the surface 
at point P are 

hR/hx Sr /S r 1 ,Sr/S0 

, y 

[grad R| I grad R] r |grad R| 

The velocity in the normal direction 
is thus 


z 



Sketch (h) 


I grad R| 


^ Sr /S r 

[grad R| 


SR/Se ^ 
Sr r2|gradR| Se 


For a quasi -cylindrical body the equation of the body surface is 


so that 


R(x,r,e) = r 



n=o 


2n0 







Since 


Scpc 

Sx 


V, 


Sr 


0 



Se 



the third term in the equation for Vn is negligible compared to the 
first two terms so that 
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Combining equations (l) and (5) yields the boundary condition for zero 
normal velocity at tbe body surface 


Sr 


^^i 
Sr Sr 


S(f>B 


I 


g 2 n(x)cos 2n9; 


n=o 


r = 1, 0 < .0 < 2ir 

(6) 


For no flow normal to the wing surface , we obtain 


ScPq/Sz Sz-u 
V 


In terms of potential functions this boundary condition becomes 


Sz 


1 ^ 
^ Se 




= V 


Sx 


r>l; e = 0,rt (7) 


Equation (7) is fiilfilled by expanding cpj_ and cpg in cosine series of 
even multiples of 0 . Equation (6) is fulfilled by satisfying the 
following equations for each term of the Fourier series as in refer- 
ence 4. 


Scp^ 


Sr 


n=o 


f 2 n(x)cos 2n0 


Sq>B 

Sr 



g 2 n(^)=°s 2n0 


n=o 


It is to be remembered that cp^ is the interference potential defined 
so that the combination possesses a cylindrical fuselage. It changes 
if the definition of the wing alone changes. All effects of distorting 
the circular cylinder are contained in cpg. 

Pressure Field Acting on Wing-Body Combination 


The first step in the present method for evaluating the drag of a 
wing -body combination is to determine the pressure distribution acting 
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on the combination. This can he done for the body very simply by means 
of the W 2 n(x) functions tabulated in reference 5* The pressure field 
off the body has been obtained in references 4 and 5 "by distributing 
multipoles along the body axis. However, the determination of the pres- 
sure field off the body can be reduced to the same amount of work as 
determining that on the body by generalizing the W 2 n(x) functions to 
W 2 n(x,r) functions. This new procedure besides vastly simplifying the 
calculation of the pressures off the body also avoids the necessity of 
integrating across logarithmic singularities involved in the original 
method. The generalization of the W 2 n(x) functions to W 2 n(x,r) func- 
tions and their use in calculating the pressure field off the body are 
described in Appendix A. The W 2 n(x,r) functions are plotted for 
n = 0 and 1 in figure 1. 

The pressure coefficients of the combination are the sums of the 
coefficients due to the body alone, wing alone, and interference. Accord 
ing to Appendix A the pressure field due to the body alone (distorted) is 
given on the body and the wing as 


^B(d) 

(tm/cr) 


I £ [s2n(==)- g 


COS 2n6 (8) 


n=o 


%(d) 

(tm/cr) 



~ ggn(x- pr+ p) 


1 

P 



x-p(r-i) 



r +1 - 



(9) 


The pressure field due to the wing alone can be determined as a pure wing 
problem in linear theory. Let he the component of this field act- 

ing on the wing and I^(y) the component acting on the body. Let 
be composed of a component; due to infinite aspect ratio plus a component; 
due to the difference between infinite and the actual aspect ratio, that 
is, a component due to the tip. Performing a Fourier analysis of the 
field around the body section and indicating the part due to the 
wing tip with an asterisk, we obtain 


^B(W) 

(tm/cp) 


00 


P2n(x) + Pan*(x) 


n=o 


cos 2n0 


(10) 
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The factor of 2 has heen introduced purely for later convenience, and the 
values of P 2 n(^) Pan*(^) represent only half the true Fourier 

amplitudes. The quantity P 2 n*(x) can he thought of as due to the image 
system of line sources and sinks representing the tip when the wing pres- 
sure field is constructed of line pressure sources in accordance with 
reference 7* 

The pressure due to the interference field is obtained in the same 
manner as in references If and 5, or analogously as that due to body dis- 
tortion. Let the normal velocity induced at the body surface by the 
wing-alone field be Fourier analyzed around the body in a manner similar 
to the wing-alone pressure field so that 


l£ 

V 



^2n(^^ + ^*2n 



cos 2n0 


( 11 ) 


Equation (ll) gives the virtual distortion of a cylindrical body neces- 
sary to conform to the wing-alone flow field. Since the cylindrical body 
is, in fact, not distorted by the wing flow field, it sets up an inter- 
ference pressure field equal and opposite to that associated with the 
virtual distortion. This opposite effect also explains the use of a 
minus sign in equation (ll) . Analogous to equations (8) and (9) , we 
obtain 


^B(i) 

(tm/cr) 


2 

P 


00 



n=o 


+ f2n*(x) - 


I 


X r 


f2n (^ ) + 


o 




( 12 ) 


^W(i) ^ 2 y^/ f 2 n(x - Pr + p) (x- Pr+P) 

(tm/cr) ^ 


1 

P 



x-pr+p 


^’an(^) + 



+ 1 
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Drag of Body and Wing in Combination 


The resultant pressxire acting on the body of the combination is the 
sum of the pressure coefficients given by equations (8), (lO), and (12). 
Let the body be distorted over the interval 0 < x < 2 , Then the body in 
combination has the drag 


(”I(C) •[ h('i) ^( 1 )] 


r ^ d0 dx 
dx 


ilh) 


With r taken as unity, ve obtain 


(B/q)B(c) 

(tin/cr)^ 


■if-ri 


g 2 m(^)^°s 2m0 


00 

r 


) i S2n(x) + PPan^^) + 


m=o 


n=o 


PP2n*(x) + f2n(x) + f2n*(x) - |^g2n(0 + fan(0 + 


WU)] W 2 n(^, <i|j cos 2n0 d0 


(15) 


The total drag of the body in combination can be expressed as the sum of 
the drags of the various Fourier components since the components are 
orthogonal; that is, the pressures due to one component cause no drag by 
acting on the shapes of the other components. 


^®/^^B(C) ^ ^°/^^B(C)2n 


(tm/cr)‘ 




n=o 


(tm/^r) 


(16) 


Then for n = 0 


J So(x) |go(x) + PPo(x) + PPo*(x) + fo(x) + V(x) - 
3 [go(i) + fo(D + fo*(o]wo (^,l^d| jdx (IT) 
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ll4- 


and for n 0 


(D/q) 


B(C) 


(tm/cr)' 


'2n 




+ PPsn^^) + 3P2n*(^) + ^2n (x) + 


t'an 


2n' 


f2n*(x) - I 


sX r 


§2^(0 + + f,„*(l) 


■2n 


2 n 


Wan 


(18) 


The pressure coefficient acting on the wing in combination is given 
by the sum of Pw(d)» %(W) # ^(i) * Consider the drag of the wing 

panel Dw(w) known from wing-alone theory, and evaluate Dw(a.) and 
as follows: Consider first the pressure components 2g2n(x 3r + 3)/3 */r" 

and 2f2n(x - 3r + 3)/3 which represent Ackeret pressures or two- 
dimensional linear-theory pressures attenuating inversely as the square 
root of r along the body characteristics. Including tip effects and 
considering both sides of both wing panels, we obtain for this component 
of the wing wave drag 


(^/q)w(i-Hi) 

(tm/cr) 


8 

3 



g2n(x-3r+3) 


SP n=o " 


Vr 


f2n(x-3r+3) 


f2n*(x-3r+3) 




dSp 


(19) 


(D/q) 


(x) 

W(i+d) 


00 

3r+3) ^ 




n=o 


Sp 




•Jr 


f2n*(x-3r+3) 

n/F 


Cr 


8n(x-3r+3,r) 

8x 


dr d(x-3r+3) 


( 20 ) 
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It is convenient to integrate first out- 
ward along strips in the characteristic 
direction and then to sum the strips 
between the leading and trailing char- 
acteristics tangent to the wing plan 
form as shown in sketch (c). 


(i) 

(tm/cr)^ 




Sketch (c) 


f2n(x-Pr+p) + f2n*(x-3r+3) 


d(x-3r+3) J 


^n(x-3r+3,r) 

cr — 


Bx 


or 


(D/q)' 




( 


w(i+d) 8 


(tm/cr)‘ 


■fix 


n=o 


X2 
xi L 


S2n(x) 


L 


So 


Cr 


^q(x,r) dr 
Sx \/r 


dx 


( 22 ) 


Consider now that part of the drag due to the pressure disturbances 
lying behind the Mach cones on which they originated and represented by 
the integral forms in equations (9) and ( 13 ) . This second component of 
the wing wave drag is 


( 2 ) 

^®/‘^^W(i+d) 

(tm/cr)^ 




x- 3 (r-i) 



+ ^*2n(0 


+ 



E 
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The volume of integration is shovn in 
sketch (d). In this sketch the posi- 
tive r direction has been reversed 
for clarity. Let us reverse the 
order of integration so that in the 
x,r,| space we Integrate over a slab 
parallel to the x^r plane and then 
sum the slabs with respect to | . 
Since WgnCx^r) = 0 if x is nega- 
tive, there is obtained 


(l>/<l)y(i+d) 


("tmAr),^ 



§2n^^^ t 


^2n(^ ^ + ^2n ^ 


d| 



r+l-|,r 


Cr ^ dr d(x-Pr+B) 
ox 


(24) 


or 


(2) 

(l>/q)y(i+d) 

(tjn/cr) 


8 

W 




=2n^^^ + ^2n^^^ ^2n 


] 


n=o 


Jf. W2n(|-r + l- |,r^cr ^dr d(|-Pr + 3) (25) 

Sp 
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The entire drag of the wing in the presence of the distorted "body is 
thus 


(tm/cr) 


(W^r)^ 




+ f2nW 




^•n(x^r) dr . 

"Jr 






+ f 


2 n 


^(x) 


dx 



Wan 



r+ 1 - 



cr 


SI 


dr d(| - pr + 3) 


(26) 


Body Shape for Minimum Drag and Drag Saving 


In this section we will derive an integral equation for the body 
shape for minimum wave drag. As noted in the Introduction, the alternate 
method of Lomax and Heaslet also. yields a solution for body shape for 
minimum wave drag. The two solutions, although not mathematically iden- 
tical, are equal to the order of quasi-cylindrical theory. The solution 
by the alternate method as well as its relationship to the present method 
is described in Appendix C. 

The drag of the combination can be expressed as the sum of the wing- 
alone drag plus the drag of a number of independent Fourier components. 

We can minimize the drag of any component or combination of components 
independently of the others; for instance, we can minimize the drag of the 
second component without altering the first component which controls the 
volume. The equation for the wave drag of the combination to be used in 
the minimization is 
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(D/q)c 

(■titt/cr)* 


°° z 

" f Z I |s2n(x) + 3P2n(x) + 3P2n*(^) + ^znU) + f2n*(x)- 

n=o ° 


I f [W5)*W5) + f2n*(5>] M^„(i|i,l)dt]-ax * 


00 2 { 

[®2n(x) + f2n(x) +f‘2n*(^) f 


n=o 


^Tifx.r) dr j 
cr ■ * '* ’ f ^ dx + 

8x vr 


(tjn/cr)^ 



n=o 


g2n(x) +f2n(x) +f2n*(x) 


dx 



j cr ^ dr d(| - 3r + 3) 
/ dl 


(27) 


The range of the integrals in the terms ms increased from Xi < x < X 2 
to 0 < X < Z with the help of simple properties of W 2 n(x^r) and ri(x,r). 
The constant An h for n = 0 and 2 for n > 0. In the minimization 
of the drag as given hy the above formula, the condition is considered in 
■which the shape of the indentation for the particular Fourier harmonic 
is at its value for minimum drag and then a variation from this 
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shape is considered. The additional drag due to the variation, 5Pp, is 
zero and is given as 


S(l>/q)p 

(im/'^r) 


“ I 

3 Z y' ^®2n(x) |^2g2n(x) + PPgnCx) + PP2n*(x) + fanCx) 


n=o 


^2n “ p Z J <ix - 

00 2 jj 

^ J' S2n(x)dx ^§2n(^)^2n ^d| + 

n=o ° ° 



Sr](x,r) dr 
Sx N/ir 


8 



Sgpn(x)lx 


n=o 



t 'll- <1(1 - 3r + 5) 


(28) 


With the inversion of the following integral 


^ g2„(x)dx^ 5g2n(^)W2n( p^l^dl-^ 5g2j,(x)dx J g^^ ( |) 0^d| 


( 29 ) 
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equation (28) Becomes 


5(D/q)c 

(im/cr)*' 


CX) I 


n=o 




i -X 




'U j 


^2n(0 + 


_. , >. 00 j 

f2n^(^)JW2n(^>lJdl|dx + I J bg^^ix) 


n=o 



8q(x,r) dr 
hx \fr 



r + 1 



Cp — dr d( I - Pr + p) 

8^ 


dx 


( 30 ) 


Now 6(D/q)^ must be zero to the first order and since 
arbitrary and not zero, ve must have for each n 
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+ PP2n(x) + PP2n*(^) +^'2n(^) 



x-g 

3 


d|- 


1 

3 



^ 2 n(^ ) + ^ 2 n*(^) 




^T|(x,r) dr 

Bx ■Vr 



0 


3 r + 3 ) 


= 0 ( 31 ) 


This integral equation defining the optimum body distortions for each 
Fourier component can be written 


g2n(x) = ll2n(x) + ^ 


^ 2 n 


U)W, 


X- 


2 n 




( 32 ) 


wherein 
ti2n(^) = 


i[ 


3P2n(x) + 3P2^*(^) + ^2n^^) + ^‘2n*(x) 


itA- 





W 2 n 



-r + 


l_x 5 n(i:^iSi£)dr d(f - 3 r+e) - f^° 

^ 31 3 x ■>/? J 
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^2n(^) + f2n^(^) 



( 33 ) 
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The nonhomogeneous term h 2 n(x) must be evaluated for the particular wing 
under consideration, the values of PgnCx), P 2 n*(x), fznlx), and f 2 n*(x) 
following from a knowledge of the wing -alone flow field. Equation (32.) 
is then solved by nvanerical means to determine the best distortion of the 
body for the Fourier harmonic under consideration. The actual distortions 
are proportional to the wing thickness because of the use of quasi - 
cylindrical theory and the application of the boundary conditions on the 
r = 1 cylinder. The wing thickness must be sufficiently thin, therefore, 
for the distortions to be within the compass of quasi-cylindrical theory. 
For a wing of greater thickness an auxiliary condition can be imposed on 
the allowable magnitude of the distortion and the minimum problem resolved 
subject to such a condition in a manner similar to that of Lomax and 
Heaslet (ref. 3)* 

A simple mle for the drag saved can be determined by evaluating the 
drag of the combination at the minimum-drag condition. Under this condi- 
tion the drag is 


(tm/cr) 


00 2 ^ 

f Z / SanW -g2n(x)+| J g2n(^)Vsn 

n=o ° '^x 


I - X 


>1 d| 


dx + 


(P/q)y(y) 

(tm/cr) 



^2n(^) +^2n*(^) 



cr 


^•n(x,r) ^ 
Sx 'Jr 


00 



n=o ° 




^2n^^) f 2 n 



1 



dr d(5 - pr + p) 




The last three terms represent the drag of the wing panels with a cylin- 
drical body, in which case the drag of the panels is the total drag of 
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the winged part of the combination. The first term thus represents the 
drag saved. If the order of integration of this term is reversed, we 
obtain 



It is readily seen that the drag saving is simply the drag of the dis- 
torted body without any wing panels, and that this drag is the same 
whether the distorted body moves in one direction along its axis or in 
the other direction. If the combination is operated at a Mach number 
other than the design Mach ntanber, the drag saving is less that that 
given by equation ( 35) ♦ 

♦ 

It is possible to build up the potential fields due to the addition 
of a quasi -cylindrical body to a wing, 9j^ and cpg, by distributing sources 
and multipoles along the x aocis as well as by the W 2 n(x,r) method. 

Any set of solutions to the wave equation in cylindrical coordinates 
that is of the form of a function of x and r mtiltiplied by cos 2n6 
can in principle be used for the purpose. In reference 4 a method is 
given of relating the body distortions to the multipole strengths for the 
particular set of solutions used therein. The same method is used in 
Appendix B to show the relationship between body distortion and multipole 
strength for the multipoles of reference 3* Appendix C the equivalence 
between the quasi -cylindrical and multipole methods is used to derive, the 
body shapes of a minimum-drag wing -body combination in a form amenable to 
analytical as well as numerical calculation. 

II - CALCULATIVE EXAMPLE OF SWEPT -WING AND BODY COMBINATION 
WITH SONIC LEADING AND TRAILING EDGES 


An interesting example of the general theory is the case of a swept- 
wing and body combination with sonic leading and trailing edges. By dis- 
torting the body properly it appears possible to reduce the drag rise 
that occtirs as a swept -wing airplane approaches the tfech ntimber for vhich 
its leading edge becomes sonic and thereby to increase the Mach number 
at ■which a given swept-wing airplane loses its sweepback effectiveness. 
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The supersonic area rule (refs. 8 and 9) gives an approximate idea of 
the shape of the distortion due to the axially symmetric harmonic . 


As 

of the' 
M = -J 2 . 


an example, a minimimn-drag computation is made for t'wo harmonics 
swept wing and "body shown in sketch (e) . The calculation is for 
If a combination at some other Mach nmher is desired, it can 
f be first transformed to M = n/I” 

in accordance with Appendix D and 
then the answers should be trans- 
foimied back to the original Mach 
number. The results of the calcu- 
lation for swept wings with sonic 
edges are plotted in nondimens ional 
form applicable to any Mach ntunber. 
A biconvex airfoil section has been 
used -with a thickness distribution 
given by 



Zu = 2% 


1 - 


x- r 

c 


= tmT][x - (r-l),r] 


(36) 


on the wing and = 0 off the plan form. The first major step in the 
calculation is to evaluate the various components of h 2 n(x) as given by 
equation (33)* The relative magnitudes of these components are in direct 
proportion to their importance in prescribing the proper shape for minimum 
drag so that considerable insight into the problem is obtained in this 
step. The second step is to solve the integral equation for the minimiz- 
ing shapes, and the third step is to evaluate the drag saving, vol\ime of 
cutout, pressure- distributions, etc. 


Velocity and Pressure Amplitude Functions 


The first four terms of equation ( 33 ) relate to the pressure and 


normal velocity fields produced at the 



Sketch (f) 


surface of the body by the wing- . 
alone flow field. For the present 
example consider the ■wing alone to 
include the area blanketed by the 
body when the leading and trailing 
edges are extended as straight 
lines. The pressure and velocity 
fields of the wing alone can be 
built up by superimposing line 
pressure sources of the type given 
by R. T. Jones in reference 7* The 
u, V, and w fields of such a line 
source, as pictured in sketch (f). 
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are given in reference 10 for 


M = as 


u 

V 


-R.P. 


me 




cosh 


-1 




+z* 


(37) 



- = R.P. - cos"^ 
V ^ 


yy' +z^ 


Vy®+z2 


(39) 


vherein 


X* = X - myj y» = y - mx; 


= zVl-m^ 


For a sonic leading edge, let m approactL unity and ottain the limits 
of equations (37)^ ( 38 ), and (39) as 


u __ G n/x^ 

’ ’ ' * (x-y) 


(1*0) 


V £ 

Y ^ It 


(x-y) 


cosh"^ 


*7y2+z2 


(i^l) 


w 

V 


€ -1 
— COS 

Jt 




y(y ^x) + z2 

(x - y)Vy^+z2 




(»^2) 


From these equations the normal velocity and pressure distribution acting 
on the cylinder y2 + = 1 can easily be determined. The velocity 

normal to the cylinder was computed for 15 ° intervals around the cylinder 
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and then Foxirier analyzed using numerical methods for the pair of line 
sources 1-1 shown in sketch (g) as well as for a pair of tip sources 2-2. 
The coefficients so calc\ilated correspond to the Fourier series 


— = fog(x) + f2g(x)cos 26 + . 

Ve 


( 1 ^ 3 ) 


= f0c;*(x) + f2o*(x)cOS 2e + . . . 


Ve 




Mach cone intersections 

with body 


wherein the asterisk refers to the tip sources. The values of the 
Fourier coefficients are shown in figure 2. The constancy of the fo^(x) 

coefficient is noteworthy since 
the normal velocity distribu- 
tion varies with x. It is 
apparent that the tip has a 
minor influence . 


To obtain the velocity - 
amplitude functions for the 
wing as a whole as defined by 
equation (ll), it is necessary 
to superimpose the effects of 
the sources and sinks that form 
the wing in accordance with 
sketch (h) . The effect of the 
infinite number of small sinks 
between the leading and trail- 
ing edges is summed by inte- 
gration. The resultant values 
of fo(x), f2(x), fo*(x), Po(x), 
p^Cx), and Pq*(x) are given in 
figure 3. 

The main terms of equa- 
tion (33) 3-^^ readily deter- 
mined. The final term involves 
numerical integration using the 
values of the W2n(^^^) 
tions in figure 1. The sixth 
term can he ascertained in 
closed form. Since from equa- 
tion (36) 



Sketch (g) 



Sketch (h) 


r\{x - r +l,r) = 2 


x_-r 

c 


1 - 


x-r 

c 
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we obtain 


^Ti(x,r) ^ 2 


2(x-l) 


(46) 


Thus 


o 

/ 6q(x,r) dr , 

2(x-l)' 

Si 6x -Jr 

L J 


i^-1) (47) 


The term involving a double integral can be reduced to a single integral 
with the aid of equation (46) 


J' J' W 2 n(^ - r + 1 - x,r) ^ |^ (^ - r + l,r)dr d(| - r + l) 

Sp 


ff W 2 n(l - r + 1 - x,r) 

2(1 -r)" 

bp 

L ^ J 


dr dU - r + 1) 


,2 


d(| - r + 1 ) j W 2 n(^ - r + 1 - x,r)dr 


= 2 


Wsnd 


+ 1 - x,r)dr 


(46) 


With the following definition 


r2n(x,s) = J W2n(x,r)dr 


(49) 
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we achieve the desired result 


f f Wsnd - r + 1 - x,r) c dr d(| - r + l) 

^ ^Sp ot 


= (l - -x + l,s)d5 


(50) 


Solutions of Integral Equations for First IVo Harmonics 


The solutions of the integral equations represented by equation (32) 
have been obtained both by iteration and by simultaneous equations. If 

(o) 

a first approximation to gg^Cx) is known, say ggn(x) , then a second 
approximation can be determined from the integral equation as 


(1) I (o) 

g2n(x) = hgnCx) J S2n(^)W2n|x-l jd^ 


(51) 


The process is repeated until it converges (or diverges) . The conver- 
gence is greatly aided by a close first approximation. If a means of 
solving many simultaneous algebraic equations is at hand, the integral 
equation can be expressed as a number of simultaneous equations by writ- 
ing the integral approximately as a summation. The more equations, the 
more accurate the answer. The rate of convergence is decreased by 
increasing the interval 0 < x < Z or by increasing the order of the 
harmonic . 

For the n = o harmonic the solution was started by using ho(x) as 
a first approximation. The component terms of hod) are shown in fig- 
ure 4. It is clear that the terms fo(x) and Pq(x) are equal and oppo- 
site and the tip effects are small, so that the effect of the body dis- 
tortions sending waves directly along the Mach cones is the dominant 
effect. This fact leads subsequently to a simple approximate theory. 

In the first calculation, the body was distorted only over the interval 
of the root chord. The subsequent approximations to go(x) were obtained 
by iteration in accordance with the equation 


go^'^''’^^(x) = ho(x) 




c+i (0) 

go(^)Wo|x-l|d| 


(52) 
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The successive approximations converge so rapidly as shown hy figure 5 
that the initial approximation could have heen used as the basis of a 
good approximate solution. In the second calculation, the possibility 
of distorting the body in front of the wing was considered with the 
result that the calculated distortion here is very close to zero, as 
shown by figure 6. For a general wing this corresponds to no distortion 
in front of the Mach rhombus enclosing the wing, a condition that follows 
readily from a consideration of momentum transfer through a large volume 
enclosing the combination. The effect of distorting the body downstream 
of the wing was investigated with the result that some distortion in this 
region is effective in further reducing the drag. However, the additional 
gain calculated for the present example is small. Thus, the significant 
interval for the purpose of the present example is the root chord, and no 
large change in the shape for minimum drag occurs over the interval of 
the root chord by considering a distortion interval greater than the root 
chord . 

A solution for the integral equation for n = 1 was attempted in the 
same way as that for n = 0, but it was soon apparent that the rate of 
convergence was relatively slow. For this reason a better initial approxi- 
mation was sought by solving a system of 11 simultaneous algebraic equa- 
tions using the trapezoidal rule to evaluate the integral. The approxi- 
mation so obtained is shown as ^(x) in figure 7. Several iterations 
then produced convergence . 


Body Shape, Drag Saving, and Pressure Distributions 


The solutions to the integral equations, go(x) and g^{x) , give the 
shape of the body for two harmonics from the integral of equation (l). 




goU)<il + (tm/c)cos 20 



(53) 


The shapes of the indentations for the n = 0 and n = 1 harmonics at the 
wing-hody junctures are shown in figure 8 for unit wing thickness. It 
is to he noted that the n = 0 indentation is similar to the wing profile. 
However, because of the finite chord length, the indentation does not 
come up to the full body radius at the wing trailing edge, a fact in 
accordance with the supersonic area rule. This point will be further 
elucidated when we consider effects of plan form on the n = 0 indenta- 
tion.. The n = 1 indentation is much larger than that for n = 0, an 
effect that tends to discourage the use of the higher harmonics. For 
the 5-percent-thick wing postulated, the indentation due to the second 
harmonic is so large that the assumption of a quasi -cylindrical body is 
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being stretched. The body cross section given by equation (53) incor- 
porating both harmonics is shown in figure 9. The change in shape behind 
the trailing edge does not contribute much to the drag reduction as will 
be discussed later, and we therefore have the alternatives of incorporat- 
ing the n = 0 afterbody, n = 1 afterbody, the afterbody using two har- 
monics, or a cylindrical afterbody as in figure 9. The percentage of the 
wing volume removed up to the trailing edge is 69 percent and up to the 
end of the Mach rhombus (x = 8) about 101 percent. Body volume is 
removed only by the n = 0 harmonic . 

The slope of the indentation at the wing leading edge can be deter- 
mined from the integral equation with the aid of the fact that 
must be nearly zero in front of the Mach rhombus (see Appendix C) . The 
discontinuities in ggnCx) must coincide with the known discontinuities 
of h2n(x) at X = 1 because the integral of equation (32) is continuous. 
Hence , 


ggn^^^) “ tL2n(l'*’) - ti2n(l") 


( 5 ^^) 


The only term contributing to the discontinuity is the next to the last 
term of equation (33). With the help of equation (^7)> we obtain 




Panels Body 
I 767 1 0 ~ 
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1 ) 


‘n 


(55) 



Panels Body Harmonics 


A + B 
A+B+C 


5.48 

0 

n = 0 

5.67 

0 

n = 0J 


(2) Wing-cylindrlcal-body combination 



Panels Harmonics and 


A+B + C+D 

A+B+C+D+ 

E 


Panels 

Body 

end 

body 

intervals of 
distortion 

3.60 

79 

4.39 

n=0,l<j^<8 

-.49 

2.61 

2.12 

n=0, l<-|-<8 
pa . 

n=l,l<j3<4 


(3) Wing -distorted -body combination 
Sketch (i) 


Since each harmonic contributes 
a finite amount, the indenta- 
tion for all the harmonics is 
initially radial in the 0=0 
and 0 = It planes. 

The net drag saving due 
to body distortion has been 
computed by evaluating the 
drag of the distorted body 
alone in accordance with equa- 
tion (35) and also from the 
pressure distributions acting 
on wing and body in combina- 
tion. The distributions of the 
drag along the wing panels due 
to two harmonics of interfer- 
ence and distortion, shown in 
figure 10, have been integrated 
to obtain the net panel drags. 


The drags of the panels and body in combination for various conditions are 
summarized in sketch (i) which is to be studied in conjunction with 
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figure 10. Part (l) of the sketch shows the drag with the panels as part 
of the wing alone. Part (2) of the sketch illustrates how the addition 
of the cylindrical body has the effect of reducing the panel. drag in the 
proportion to 5*^7. (it has been ass-umed that the interference drag 

for harmonics higher than the n = 1 harmonic is negligible.) Part (3) 
of the sketch shows the drag saving due to two harmonics of distortion. 
Incorporating the n = 0 distortion into the body reduces the drag of the 
combination by 22 percent (of ^.6j) and the addition of the n = 1 harmonic 

has the effect of reducing the drag of the combination another k2 percent 

for a total of 6h percent. It is to be noted that distorting the body 

has the effect of reducing the drag of the wing in the presence of the 

body while increasing the drag of the body in the presence of the wing. 

The wing has negative wave drag when two harmonics of distortion are 
incorporated into the body. The net gain in drag is the result of reduc- 
ing the wing-panel drag about twice as much as the drag of the body is 
increased. It is interesting that the section drag coefficients given by 
figure 10 for the harmonics of interference and distortion attenuates 
along the characteristics approximately inversely as the square root of 
the radius, so that drag savings can be realized for aspect ratios greater 
than that of the present example . 

On the body the drag loading can be defined from equations (I7) 
and (18) as 


dx\qSHy^(C) 

(Wc)^ 


n ^ ^ ^ f* 

+^’2n(^) +^2n*(^) " 

n=o 



S2n(^) + ^2n (0 + 


W2n(x- l)d|| 


( 56 ) 


This quantity, plotted in figure 11, for the n = 0 and n = 1 harmonics, 
shows the distribution of drag along the body. This drag results from 
the pressure field arising from the body distortion and by distorting the 
body in an existing pressiire field. It is noteworthy that extending the 
interval of distortion aft of the wing trailing edge causes only minor 
reductions in wave drag for the first harmonic. An approximate calcula- 
tion shows this also to be true for the second harmonic . 

The pressure distributions acting on wing and body both in combina- 
tion can be readily determined, the distribution in the wing-body juncture 
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"being given by 


P 

(‘tm/c) 


n=r 




f2n(0 + f2n*(t) 


W 2 n(x- l)dl| 


S2n(x) 


J g2n(^)W2n(x-0d| 


P2n(^) 


+ 


P2n*( 



( 57 ) 


The first component represents the effect of interference between wing 
panels and cylindrical body, the second component the effect of body 
distortion, and the third component the effect of the wing-alone field. 

The pressure distri"bution at the wing-body juncture for a cylindrical 
body is shown in figure 12(a) . The wing-alone component of this pressure 
distribution has an infinity at the leading edge of the wing -body juncture. 
Because the body is a perfect reflection plane as far as this point is 
concerned, the pressiire coefficient is the same as that at the apex of the 
wing alone, namely, 8/it. The finite contributions of a finite number of 
harmonics of interference cannot overcome this infinity. However, we can 
easily establish the curve for the combination near the leading edge by 
fairing it into the known value. It is to be noted that the resultant 
pressure distribution has a favorable pressure gradient for the boundary 
layer. In figure 12(b) the pressure distribution for a distorted body 
has been determined by adding the pressures due to two harmonics of dis- 
tortion to those for the cylindrical body. The resulting distribution 
has an unfavorable gradient. The effects at the wing -body juncture are 
typical of the wing as a whole since the pressure distributions on the 
span are essentially those of the root chord with an inverse variation 
as the square root of the radius (fig. I 3 ) • They are also typical of the 
body since they are simply multiplied by a cos 2n0 factor. For a 
5-percent-thick wing the complete -pressure field at the wing-body juncture 
shows a pres sure -coefficient rise from about -0.20 at the leading edge to 
about 0.25 at the trailing edge. This rise of Ap/q = 0.45 (Ap/p = O. 63 ) 
is less than the sudden pressiire rise, Ap/p, of 0.9 necessary to separate 
a turbulent boundary layer as given by Bogdonoff and Kepler (ref . 11) . 

It is thus probable that turbulent separation will not occur. 0 


Approximate Solutions for Axially Symmetric Harmonic; 
Effect of Plan Form 


On the basis of the relative magnitude of various terms of ho(x) in 
the integral equation for go(x), equation ( 32 )> it is possible to obtain 
simple approximate solutions for the indentation shape and drag saving 
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for the n = 0 harmonic. The relative magnitude of the various terms 
depends primarily on the wing chord, or more specifically on c/ 3 r, so 
that a short-chord theory and a’ medium-chord theory can he distinguished. 
In the analysis we will confine the indentation to the root chord and 
let 3 Ae be unity or greater so that no tip effects occ\ir. The Mach 
nianber need not be although the analysis is carried out for this 

Mach number. Under these conditions the 'effects of interference and the 
wing-alone pressure field represented by fo(x) and Pq(x) are equal and 
opposite for this particular case, and equations ( 32 ) and ( 33 ) can be 
rewritten 


go(x) 

= ho(x) + 1 y' go(^)WoU-^ l<3-5 

(58) 

ho(x) 

= boj^(x) + ho2(x) + hog(x) 

( 59 ) 

with 



boi(x) 

^ r^°c ^ 

Jg^ sfr 

(60) 

bogCx) 

bp 

(61) 

CO 

"if 

(62) 


o 


The components of ho(x) involving the W 2 n(x,r) functions are concerned 
with departtire of the pressures due to body indentations from the values 
that would be calculated by linear two-dimensional theory, and as such 
these components assume importance for medi\mi to long chord lengths only. 

Short-chord theory .- By neglecting ho 2 (x) and hQg(x) compared to 
boi(x) , we obtain an approximate simple theory that represents the solu- 
tion for vanishing chord length. Under this assumption the pressures 
acting on both body and wing are simply the linear-theory two-dimensional 
pressures based on the body distortion attenuated inversely as the square 
root of the radius, This simple physical model contains the primary 
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features of the flow due to the n = 0 harmonic, and the medium-chord 
approximation can he considered a slight modification of the simple model 
to account for three-dimensional effects on the pressure field. On the 
basis of the short-chord approximation, equations (58) to (62) yield 

go(x) =hc(x) = - i <= (63) 


utilizing the result of equation (47), we obtain 



Thus the shape of the indentation is given as 

(l-rg) = ^ ^ 1)ti(x,1) (66) 


The shape of the indentation is similar to the wing section^ hut the 
volume of the indentation is usually less than the volume of the wing. 
In fact the ratio of indentation volume to wing volume is 


Vi nTs - 1 

Vy S - 1 


(67) 


The volume ratio decreases as the span increases although the actual vol- 
\ime removed increases . 



NACA RM A55B07 


35 


The drag saving calcvilated using equation (35) is 



(68) 


Expressed as a drag coefficient based, on the area of the wing panels, 
the drag saving is 


^ _ 32 (s^-l) 

(tjn/c) (‘fcm/c) (‘v/s + 1) 


The coefficient of drag saving thus increases as the span increases but 
so does the drag of the wing alone. The effect of chord length will sub- 
sequently be discussed. 

Medium-chord theory .- It was possible by neglecting three-dimensional 
effects to obtain very simple results for the first harmonic . However, 
for longer chord length than those to -vdiich the simple theory applies > 
account must be taken of the terms containing W2n(x,r) functions. The 
value of hoi(x) remains as before. In the evaluation of hogCx) a simple 
assumption permits an approximate determination to be made readily. This 
assumption, expressed in physical terms, states that the three-dimensional 
effects represented by the Wo(x,r) function themselves attenuate along 
the body characteristic in accordance with a simple approximation similar 
to that for the two-dimensional effects. Mathematically the approximation 
can be expressed as 


Wo(x,r) 


«« Wo(x,l) 


Wo(0,r) 

Wo(0,l) 


(70) 


The validity of the assumption can be measured by how closely the quantity 
Wo(x,r)WQ(0,l)/Wo(0,r) is invariant to changes in r. A plot of this 
quantity for various values of r is given in figure l4 to show that the 
approximation is a good one. The value of W2n(0,r) obtained from Laplace 
transform theory is 


= + 3 _ l6n^ - 1 

8n/F 8r®^^ 


W2n(0,r) 


(71) 
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so that 


Wo(x,r)« (^ + _ 


(72) 


From equation (72) we obtain the result 


(X) , j- f r 3 1 (, 3 , 

’^^O ^O "^1 8Wn(0,l) \r^/2 1 . 3/2 cy 


IO2 


Integration yields 




(*«/s - 1 ) ( 3 ^/T + 1 ) 


2 (x-l) 


Y(c-x + 1)-^Z (c-x + l) 




1 < X < c + 1 (7k) 

with the following simple numerical quantities given in figure I 5 : 


Y(x) = / Wo(l)d|; r = 1 


(75) 


>x 


Z(x) = / IWo(l)dl; r=l 


( 76 ) 


Likewise the value hog(x) with the help of 


fo(l) = - I 


2(1 - 1 ) 


; 1 < I < c + 1 


(77) 


can easily be expressed in terms of the foregoing numerical quantities 
as 


bos(x) 



2 (x - 1 ) 

c 


Y (x-1) + I Z (x-1) 


} 


(78) 


With a good approximation to ho(x) including three-dimensional 
effects, we can now solve the integral equation for the body shape. In 



MCA RM A55B07 


37 


evaluating the integral we will approximate go(x) by bo(x) and apply 
the trapezoidal rule to the interval 1< x <c + 1 broken into n parts 
of length c/n. Then 

n 

go(xi) = ho(xi) + -^ ^^J^lo(^j)Wolxi-lj| (79) 

j=o 


vhere 






0 < j < n J 


(80) 


The example already solved exactly for the interval l<x<c+lis used 
in figure l6 to measure the accuracy of the approximate methods. For a 
chord-radius ratio of 3 the short-chord theory is already inaccurate but 
the medium-chord theory is in good agreement with the exact solution. 

Some improvement in accuracy can be obtained by iterating equation (79), 
but for a more exact solution a better integration formula is needed - 
one recognizing the slope discontinuity in Wo|xj^-li| at xi = |i. The 
drag saving evaluated by the trapezoidal rule is 

S J=0 


Effect of chord-radius and semi span -radius ratios .- The effect of 
the geometric variables on the body indentation, drag saving, and volume 
of indentation have been calculated using the medium-chord theory to 
study the limitations of the short -chord theory. The drag saving is 
shown in figure 17 ns a function of s/a and c/pa. The saving increases 
as the wing si)an increases and increases slightly as the chord increases. 
The short-chord theory gives a good rough measure of the drag saving up 
to values of c/|3a of the order of U. The volume of the indentation as 
a fraction of the wing volume is given in figure l8. It is seen that 
short-chord theory gives accurate volumes only for very short chords. 

Some actual indentation shapes are given in figure 19. The main feature 
of the shapes is that for the longer chord lengths the indentations do 
not return to the original radius. This result would be predicted by the 
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supersonic area rule. Those that do tend to return do not have nearly 
as much volume as the wing panels. Thus, as the design Mach number is 
increased from unity, the leading edge remaining sonic, the body volume 
to be removed decreases from the wing volume at M = 1 to a fraction of 
the wing volimie at higher supersonic speeds. If distortion is made over 
the entire length of the Mach rhombus, the transonic area rule gives an 
indentation volume equal to the volume of the wing panels at M = 1 and 
the present method gives 101 percent of that voltime at M = n/^. 


CONCLUDING REMARKS 


The Nielsen-Pitts theory of wing -body interference has been applied 
to the problem of minimizing the wave drag of wing-body combinations by 
utilizing radial body distortion. The method, applicable to combinations 
with quasi -cylindrical bodies and horizontal planes of symmetry, deals 
directly with press\ire distributions and body shapes. An integral equa- 
tion is derived for the body shape for minimum drag for each Fourier 
component of the radial body distortion. To demonstrate the use of the 
method, it is applied to the wave-drag reduction of a 45° swept-wing and 
body combination utilizing a wing chord of three body radii and a span 
of six body radii. The leading and trailing edges are taken as sonic, 
and the wing section is biconvex. The numerical results for this example 
support the following statements: 

1. The axially symmetric harmonic of body distortion causes a 
reduction in wave drag of about 22 percent of the drag of the wing panels 
on the undistorted body while the second harmonic causes a further reduc- 
tion of 42 percent. 

2. The over -all wave-drag reduction is the result of reducing the 
wing-panel drag by about twice as much as the body drag is increased. 

In fact, with two harmonics the wave drag of the panels is approximately 
zero . 


3 . The axially symmetric harmonic yields results qualitatively in 
accordance with the supersonic area rule in that the body indentation 
does not return to the full body radius at the trailing edge of the wing- 
body juncture. The volume of indentation is only about two thirds of 
the wing volimie for an interval of distortion of the root chord. 

4. The maximum distortion of the cross section at the wing -body 
juncture due to the second harmonic is about five times that due to the 
axially symmetric harmonic although no volxame is removed by the second 
harmonic . This limits the wing thickness for which the second harmonic 
can be fully utilized to about 5 percent or less. 
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5 . The interval over ■which the body distortion ■was effective in 
producing reduction in ■wave drag is that of the root chord although some 
fvirther small saving could be obtained by distorting the afterbody. 

6. The pressure distributions produced on the wing by body distor- 
tion, although favorable for the wave drag, are u nf avorable for the 
boundary layer . 

7 . A very simple short -chord theory can be obtained for the axially 
symmetric harmonic by neglecting three-dimensional body effects. An 
approximate treatment of these effects yields a simple medium-chord theory. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Feb. 7 , 1955 
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APPENDIX A 

PRESSURE FIELD DUE TO DISTORTED BODY AT ANY SUPERSONIC SPEED 


Let us determine the pressure field due to a "body distorted in 
accordance with equation (l) 


^r 1 


hx J 


r=i 



00 

^ g 2 n(^)‘=°s 2n0. 
n=o 


(Al) 


We will consider M = for the moment. The potential cf>g is to 
satisfy equation (Al) on the body and to fulfill the wave equation. By 
analogy with equations (12) and (l6) of reference h, 


^B 


V 



00 



n=o 


G2n(s)K2n(sr) 
®^2n* ) 


cos 2n0 


(A2) 


wherein 


^snC^) ~ ^[g2n(^)l (A3) 

and K 2 n(s) is the modified Bessel function of the second kind. For the 
case r = 1 it was possible in reference 4 to express the inversion of 
equation (A2) in terms of characteristic functions independent of the 
shape of the bumps of the body. Thus in terms of the characteristic 
functions 


W2n(x) H L‘^1 


^2n ( ^ ) 
(_K 2 n* (s) 


+ 1 


(Alf) 


it was shown that on the body 


p ^ _ 2(^tp/Sx) ^ ^ 

V 


®2n(x) - y' g2n(^)^2n(x- l)d| 


cos 2n0 


n=o 


(A5) 



NACA RM A55B07 


4l 


It is possible to determine the pressure at any point off the body 
in a manner similar to that for points on the body. This is accomplished 
by generalizing the characteristic fiinctions to functions of two variables 


W2n(x,r) = L 


= T-i 


s(r-i) ^ 2 n(sr) 
^2n (s) 


(a6) 


In Laplace transform notation the pressure coefficient is 


= -2 (i) J Oan(s) cos 


n=o 


K2n’(s) 


00 

= -2 (I) e-=<'-^>G^„(s) 


n=o 


i(r-i) Kgn(sr) ^ ^ 

Kan'Cs) ^ 'Jr j 


cos 2n0 (at) 


With the aid of the following equalities from Laplace transform theory 


. -1 


-s(r- 


^2n(s) J - 


g 2 n(x - r + 1) 


(A8) 


-1 


K2n(sr) ^ _e 


-s(r-i) 


Kan'Cs) 


•Tr* 


= W 2 n(x - r + l,r) 


(A9) 


L GgnCs) 


K; 


■ 2 n 


(sr) 


K 2 n’(s) s/r" 


px-r+1 

= J S2n(i)W2n(x- r + 1- 5,r)d| (AlO) 


the final result for the pressure field on or off the hody is 


n=o 


X - r + 1 - l^r)d| 


cos 2n6 

(All) 
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For any Mach number the pressure coefficient is 

g^„(x-|3r.|3) - 




n=o 


n/F 


px-pr+p / t \ 

- I f g2n(^)Wan(|-r+l-|^^/^ 


cos 2n0 


In this equation the functions g 2 n(^) taken to be the actual ampli- 
tudes of the body distortion and do not depend on Ifech mmiber. If the 
combination is distorted as in the Prandtl-Glauert transformation (see 
Appendix D) , then the g 2 n(x) functions will change with Ifech number. 

To make computations using equation (A12) it is necessary to evaluate 
the integral numerically using tabulated or plotted values of W 2 n(^^^) • 
Plots of Wo(x,r) and W 2 (x,r) sufficient for the purposes of this report 
are included in figure 1.. 
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APPENDIX B 

RELATIONSHIP BETWEEN SHAPE OF QUASI -CYLINDRICAL BODY 
AND STRENGTHS OP AXIAL MULTIPOLE DISTRIBUTIONS 


The method used in this report deals exclusively with the flow field 
external to the wing-hody comhination and requires no continuation of the 
solution inside the body. However, it is intuitively clear that any 
quasi -cylindrical flow can he huilt up hy superimposing distributions of 
sources and multipoles along the axis of the body. For the particular 
set of sources and multipoles used in reference 4 a set of characteristic 
functions, the M2n(x) functions, were defined and tabulated to allow the 
axial distributions of multipole strengths to be determined for a given 
body shape. (The problem of finding the body shape from the multipole 
distribution is a direct, although laborious, mathematical process.) The 
same characteristic functions are applicable to the particular set of 
multipoles used in reference 3. These multipoles written in several 
forms as follows 


cos 2 n 6 cosh (2n cosh"^ p-) 
n/x^ - r^ 


cos 2n8 
a-Jx^-r^ 



= cos 2neL~^[K2n(rs) ] 


(Bl) 


are distributed along the axis in strength A2n(x) . The multipoles are 
started along the body axis at x = 0, as shown in sketch (j) so that 


r,A2„(x) 



Sketch (j) 
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the body distortion starts at x = 1 for M = The multipoles are 

not started at negative values of x since A 2 n(x) would not then have 
a Laplace transform when x < 0. The Laplace transform of the potential 
due to the multipoles is 


9M(s,r,0) 


= - ^ ) A 2 n(s)K 2 n(sr)cos 2n0 


(B2) 


n=o 


and the potential in the physical space is 


i?M(x,r,0) 


00 



n=o 


(cos 2n0) 





cosh ( 2n cosh"^ — ^ 


7(x- 0^- 


d5 


(B3) 


Let Cf^ be equal to cpg produced by a distortion of the body given by 
equation (Al) . Then, according to equation (A2) we have 


TM(x,r,0) = cp^(x,r,0) =V(— jL 


''B 


"tm 


cr 


G2n(s)K2^(sr) 


2n' 


n=o 


^^2n * ( s ) 


cos 2n0 


(b4) 


Comparison of equations (B2) and (b 4) establishes the relationship 
between body shape and multipole strengths in the Laplace space as 


G'2n(®) ~ " 2^211* ( s)A2n(s) 

2«V(tm/cr) 


(B5) 
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Form the Faltung integral corresponding to equation (B5) with the help 
of the relation 


L 


-1 



= L"^ [K2n(s)] 


cosh(2n cosh“^x)^ 

n/x2 - 1 


X > 1 


= 0; X < 1 


Thus 


§2n(^) ~ 


1 

2JtV(tm/cr) 



cosh(2n cosh'^x- I) 
■Jix-if-i 


dl 


1 

2rtV(tm/cr) 


J Asn (t - l) 


cosh(2n cosh“^ x - | + l) 
J(x - I + l)^ - 1 


dg 


(b6) 


Equation (b 6) gives a method of determining the shape of a quasi - 
cylinder corresponding to a given distribution of multipoles along the 
axis. The reverse problem of determining the multipole distributions 
for a given body shape is not as simple. This problem will arise when- 
ever a minimum-drag wing-body combination is operated away from the 
design point. Equation (B5) is written in a form to facilitate the solu 
tion of the reverse problem 


e'^ Asnis)--2iN(^^ 


G2n(s) 

se^Kan’Cs) 


(BT) 


The occurrence of K 2 n'(s) in the denominator makes the determination of 
the inverse transform difficult if at all possible in terms of known 
functions since the zeros of K 2 n*(s) in the complex plane are generally 
unknown. However, if we define the following function by its Laplace 
transform. 


^2n(^) — ^ 


L se K 211 (s) 


(b8) 
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then the solution to the reverse prohlem can be written down 


A2n(x - l) = -2itV 


(i)f 


g2n(0M2n(x- l)d| 


(B9) 


Tabulated values of Mo(x), MgCx), M 4 (x), and MeCx) for 0 < x < 4 are 
included in reference 5 for nmerical work. Equation (B9) thus allows a 
simple calculation of the multipole strengths for a given shape. The 
functions M 2 n(x:) have logarithmic singularities at the origin. Thus 
whenever S 2 n(^) finite, so is A 2 n(^"l)> and if g 2 n(^) ^ singu- 

larity, so will A 2 n(x - l) . 
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APPENDIX C 

ALTERNATE METHOD FOR DETERMINING BODY SHAPE 
FOR MINIMUM-DRAG COMBINATION 


In this appendix we will consider the solution for the hody shape 
for minimum drag hy the method of reference 3^ and its relationship to 
the present method. In reference 3^ the potential function for a general 
class of minimum-drag combinations is found without recourse to solving 
a boundary -value problem as already described in the Introduction. If 
the stream tubes corresponding to the potential are known, any stream tube 
can be replaced by a body to form a minimum-drag wing -body combination. 

To carry out the precise determination of the body shape requires the 
simultaneous solution of a pair of nonlinear differential equations. 
Practically, the shape must be calculated by some approximate method. 
However, the potential function for the minimum-drag combination is accu- 
rate to the order of linear theory. In the quasi -cylindrical theory it 
is assumed that the flow is known for a cylindrical body and wing combina- 
tion by the use of the interference methods of references Ij-, 5^ and 6. 

Then the shape perturbation from the circular cylinder necessary to pro- 
duce minimum drag is calculated harmonic by harmonic . In this method the 
boundary conditions are applied on the r = 1 cylinder. The potential 
for the minimum-drag combination will not be mathematically identical to 
that given by the multipole solution because the body boundary conditions 
in the quasi -cylindrical problem have not been applied precisely on the 
body. (For both solutions the wing boundary conditions are applied on 
the z = 0 plane instead of on the wing surface.) However, for thin 
enough wing panels the minimvim-drag body will approach a cylindrical body 
as closely as desired. Under these conditions we may therefore equate 
the potentials given by the two methods to the order of linear theory. 

In the present paper, the combination ’potential is 

9C = Tw + <Pi + C(>B (Cl) 

and in reference 3 the combination potential is 

'Pc = "Pw + ^ (C2) 

where is the potential due to the canceling multipoles. We can thus 

say that to the order of linear theory 


'i + TB = *Pm 


(C3) 
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Since from equations (B2) and (a4) the following relationships are true, 


00 

L(cPm) = - — ) cos 2ne A 2 n(s)K 2 n(sr) 

2if Z_j 


n=o 


00 

-«Pb) V cos 2 n 0 


L 

n=o 


G2n(s)K2n(si’) 

sK2n'(s) 


00 

2ne 


n=o 


F2n(s)K2n(sr) 

sK2n* (s) 


we obtain the following relationship among the variables 


^2n(®) = ~F2n(s) A2n(s) sK2n’ 

2ad/^( tm/ Cj-) 


(C4) 


The following relationships pemit us to form the Faltung integral 
corresponding to equation (C4) . 


. -1 


K2n’(s) 


= L-^ tK^nCs)] - ££SM2S.i;2^); 

7xS- 1 


X > 1 


0 


X < 1 


(C5) 


The inverse of equation (C4) is thus 


S2n(x) = -f2n(x) 


2jfV(tm/cr) 


,^°sh[2n cosh-^x - |)] 


7(x-l)^-l 


(C6) 


The first term can be thought of as the virtual distortion of the cylinder 
by the wing-alone field, and the second terra is the body distortion due 
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to the canceling multipoles. The multipole strengths are determined by 
the method of oblique planes as used in the supersonic area rule and 
described in references 3 and 8. 

The value of the shape amplitude functions, g 2 n(x), can, in prin- 
ciple, be determined analytically from equation (C6). The determination 
involves first calculating the normal velocity field at the r = 1 cylin- 
der due to the wing alone and performing a Fourier analysis to obtain 
the fan(x) functions. Then the multipole strength distributions for the 
wing alone A 2 n(l) are obtained by the method of oblique cuts which 
includes a Fourier analysis. Finally, to obtain the shape amplitude 
functions, the integrals of equation (C6) must be evaluated. 
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APPENDIX D 

PRANDTL-GLAUERT TRANSFORMATION FOR WING-BODY COMBINATION 


The calculations for the example configuration have all been made 
for a Mach number of n/ 2, but the results can be generalized to any Mach 
n\imber by means of the Prandtl-Glauert rule. Let the configuration at 
Mach number M be transformed into an equivalent M = •12 configuration 
by dividing its longitudinal dimensions by 3 and leaving its lateral 
dimension unchanged. Let the potential at corresponding points be the 
same (see fig. 20) . The two potential fields are solutions to the dif- 
ferential equations given in the figure. Since the potential is the same 
in the y and z dimensions, Bcp/Sy and Scp/Sz are unaltered so that the 
slopes of all surfaces in the streamwise direction are unaltered. Under 
these conditions the quantities at M and at M = */2 are related as 
follows : 


'M 


C y* 

^M 




®M 

= 


= 



Dm 

= A/i" 

x,y,z 

at M and 

1 

' 3 
/- 

^n(F 


4 ' 






(Dl) 


J 


dx 




r = 1 


0 = 0,jt 


(^2) 
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The calculated results at M = can be generalized to apply to 
all combinations of the family - in our case all combinations with 
sonic edges and fixed body radius and combination span - by plotting 


l-(rg/a) 


versus x/pa for body indentations, pP/(tm/c) versus x/pa 


P (tm/ c) 

with r/a as parameter for wing pressure distributions, PCp/(tm/c)^ 
versus s/a or PAe with c/pa as parameter for drag, and' Y±/Y\j versus 
s/a or PAe with c/pa as parameter for volume of cut-out. 
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(a) WQ(x,r) 


Figure 1.- Charts of characteristic functions. 
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Figure 2.- Velocity amplitude functions for source pairs. 
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Figure 3*~ Velocity and pressure amplitude functions for wing alone. 
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Figure k.- Components of ho(x). 
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Figure 5 



-Successive approximations in solution to integral equation 

for n = 0. 
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Figure 6.- Dependence of go(x) on interval of distortion. 
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Figure 7«- Successive approximations in solution to integral equation 

for n = 1. 
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Figure 8.- Contribution of axially symmetric and second harmonics to 
body indentation at wing juncture. 



Figure 9»- Shape of distorted body. 



62 


NACA RM A55B07 



A-- Panels in wing alone 

A+B Paneis on cylindrical body, n = 0 harmonic 

A+B + C Panels on cylindrical body, n = 0 and n = l harmonics 

A + B+C+D-- Ponels on distorted body, n = 0 harmonic 

A+B+C+D+E Ponels on distorted body, n=0 ond n = l harmonics 


Figure 10.- Spanwise distribution of drag on wing panels in combination 

with various bodies. 





63 



Figure 11.- Drag loading along body due to n = 0 and n = 1 harmonics of distortion. 
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( a) Cylindrical body . 


Figure 12.- Pressure distribution at ving-body juncture. 
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Figure l4.- Vatrlation of Wo(x,r) Wo(0,l)/Wo(0,r) with radius. 
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Figure 15.- Values of Y(x) and Z(x) . 
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Figure l8.- Effect of span and chord on volume of indentation for n = 0 harmonic as fraction of 

wing volume. 
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Figure 20.- Relationship between wing -body combination and equivalent combination at 





